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Abstract 

The temporal pseudoscalar meson correlation function in a QCD plasma is investi- 
gated in a range of temperatures exceeding T c and yet of experimental interest. Only 
the flavour-singlet channel is considered and the imaginary time formalism is em- 
ployed for the finite temperature calculations. The behaviour of the meson spectral 
function and of the temporal correlator is first studied in the HTL approximation, 
where one replaces the free thermal quark propagators with the HTL resummed 
ones. This procedure satisfactory describes the soft fermionic modes, but its appli- 
cation to the propagation of hard quarks is not reliable. An improved version of 
the so-called NLA scheme, which allows a better treatment of the hard fermionic 
modes, is then proposed. The impact of the improved NLA on the pseudoscalar 
temporal correlator is investigated. 
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1 Introduction 



The meson spectral functions in the different channels (scalar, pseudoscalar, 
vector, pseudovector) in the deconfined phase of QCD have been explored with 
various techniques, in particular with the Hard Thermal Loop (HTL) approx- 
imation [1,2,3,4]. 

In this paper we focus most of our attention on the pseudoscalar channel where 
we compute the spectral function in a scheme improving upon the HTL ap- 
proximation. The degrees of freedom we deal with are light quarks (massless 
for the sake of simplicity) and gluons and our results refer to the case of zero 
chemical potential (i.e. zero baryon density), a condition which is expected 
to be realized in the heavy ion experiments at RHIC and even better in the 
future experiments at LHC The calculations are performed in the imaginary 
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time formalism, taking at the end a proper analytical continuation to real 
frequencies when required. 

To introduce the subject of our investigation we first shortly revisit the free 
case, where one computes the 2-point free meson correlation function and its 
spectral density through the convolution of two free thermal quark propa- 
gators. Then we evaluate, for illustration in the pseudoscalar case only, the 
spectral function in the HTL approximation (namely using HTL resummed 
propagators for the fermionic lines in the loop), thus checking the results lately 
obtained by Karsh et al. [3,4]. 

Actually, for its phenomenological importance, much attention has been payed 
in the literature [1,2,3,4,5,6] to the vector channel. Indeed the vector spectral 
function is the basic ingredient to evaluate the dilepton production rate from 
the plasma [1] , possibly of experimental relevance. However in the vector case 
it turns out that, in addition to the quark propagator, the quark-photon ver- 
tex should also be dressed in the HTL approximation, at variance with the 
pseudoscalar case where the vertex correction vanishes. We shall later shortly 
discuss this item. 

One may worry that at the temperatures reached in heavy ion experiments 
the QCD coupling g is not so small (in fact g ~ 1) and this could cause 
problems to approximation schemes, like the HTL, based on the separation 
of different momentum scales (i.e. g 2 T, gT, T) 1 which strictly holds in the 
weak coupling regime. Nevertheless the HTL approximation (especially in its 
"improved" versions) is able to reproduce quite well the lattice data for the 
thermodynamics of the QGP phase for temperatures T > 2.5T C as reported 
for example in [7]. Hence we expect such an approximation to provide reliable 
results for other observables as well, in the same range of temperatures. 

Indeed the kinematical approximations which allow an analytical calculation 
of the quark self-energy, the basic ingredient of the HTL resummed propa- 
gator (see Appendix B), are justified only if the components of the external 
quark four-momentum are soft, i.e. k ~ gT. Hence we improve upon the HTL 
calculation of the meson correlation function by using a prescription which 
allows for a better treatment of the hard fermionic modes. In this approxima- 
tion, which was introduced in the analysis of the thermodynamics of the QGP 
[8,9,10] and is referred to as the Next to Leading Approximation (NLA), ef- 
fects arising from the hard quark modes are taken into account in an effective 
way through a correction to the asymptotic quark thermal mass. Actually we 
propose a variant of the NLA which allows for a smoother transition between 
the two regimes of soft and hard momenta. 



At each scale different aspects of the QGP physics disclose themselves. 
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Finally, as a general remark, we observe that, within the present approxima- 
tion scheme, as long as we confine ourselves to the study of the meson spectral 
functions (i.e. the imaginary part of a polarization propagator) we do not need 
to face any problem related to renormalization. 

Our paper is organized as follows. 

In Sec. 2 we recall the definitions of the basic quantities we address: the 
meson current operator, the meson spectral function and the meson correlator 
along the (imaginary) temporal direction. In Sec. 3 we display the results one 
obtains in the free case. In Sec. 4, after recalling the expression of the HTL 
quark spectral function, we use HTL resummed quark propagators to obtain 
the meson correlation function: this accounts for important effects coming from 
the interaction of the quarks with the thermal bath. Our calculations confirm 
the results previously obtained by Karsch et al. [3,4]. In Sec. 5 we examine the 
NLA scheme. Finally, in Sec. 6, we display the numerical predictions for the 
pseudo-scalar meson spectral function and for the temporal correlator G(r) 
within the two approximation schemes (HTL and NLA) and compare them 
with the free case. 



2 Thermal meson correlation function 

Adopting the same notation as in Ref. [4] we consider the following current 
operator, carrying the quantum numbers of a meson in the flavour-singlet 
channel (there is no matrix mixing the flavours): 

Jm{-iTi x) = q(-ir, x)T M q(-ir, x) , (1) 

where Y M = 1, 7 5 , 7 M , 7^7 5 for the scalar, pseudoscalar, vector and pseudovec- 
tor channels, respectively. We next define the fluctuation operator J M as 

J M (-ir, x) = J M (-ir, x) - (J M (-ir, x)} , (2) 

the average being taken on the statistical ensemble. 

The chief quantity we address is the thermal meson 2 point correlation func- 
tion: 

G M (~ir, x) = (J M (-iT, x)J ] M (Q, 0)) 

= (J M (-iT,x)jl(0,0)) - (J M (-iT,x))(jl(0,0)) 

4 £/^--""v-w-„,p), <3) 
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with r e [0,0 = 1/T] and cu n = 2mcT (n = 0, ±1, ±2 . . . )• 

It is convenient to adopt the following spectral representation for the meson 

propagator in momentum space: 

(4) 

Hence it is possible to express the thermal meson propagator in a mixed 
representation through the spectral function om- Indeed starting from the 
definition: 



J +oo 

G M (-ir,p) = - e~ lUJnT XM(iu n ,p) 

P n=-oc 



and performing the sum over the Matsubara frequencies with a standard con- 
tour integration in the complex u plane [2,11], one obtains [4]: 



n < ■ \ f +C ° a < A cosh(o;(r-/3/2)) 
GM(-.r,p) = / o &;<M",P) sinhM/2) • (5) 



sinh(o;/3/2) 



In the following, confining ourselves to light mesons (namely the pions if one 
considers the pseudoscalar flavour triplet channel), we will explore in various 
approximation schemes the behaviour of the meson spectral function. 
Indeed the survival of a meson as a bound state in the deconfined phase of 
QCD should be signaled by a delta like structure of the spectral function. 
We will see that this is not the case (at least according to the approaches 
adopted in this paper, which keep into account only the interaction of the 
fermionic modes with the thermal bath). On the other hand, for what concerns 
heavier mesons, recent lattice results [12,13,14], obtained in quenched QCD, 
do display very pronounced peaks for the charmonium spectral functions in 
the pseudoscalar and vector channels, which survive till T ~ 2.25T C . Actually, 
as recently pointed out by Shuryak [15,16], the existence of a whole sequence 
of loosely bound states above T c (not necessarily colour-singlet states) would 
entail large cross sections among the constituents of the plasma; this in turn 
would justify the adoption of hydrodynamical models for the fireball arising 
from the ultra-relativistic heavy-ion collisions; it would also explain certain 
collective phenomena, such as the elliptic flow, observed at SPS and more 
recently at RHIC. The QGP, at least for T just above T c , seems to behave more 
like a liquid with a very small viscosity rather than like a gas of quasiparticles. 
In the light of these results, which indicate the survival of mesonic resonances 
in the QGP phase, the existence of bound states in the broken-symmetry phase 
of the three-dimensional Ising model [17,18] gains considerable interest. We 
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Fig. 1. Diagram illustrating the free meson correlation function. The propagation 
lines are free thermal quark propagators. The two interaction vertices correspond 
to r matrices. 

remind the reader that, on the lattice, the 3d Ising model is related through 
a duality transformation to the Z2 gauge model [19], the center of the group 
SU(2). 



3 Free spectral functions 

We start by revisiting the computation of the free meson spectral functions for 
the different channels. We use the short-hand notation (r, x) for (t = —ir, x). 
In order to display more transparently how one gets the final result, we write 
explicitly the spinor indices. We have to compute the expression 

(J M (r, x)f M (0, 0)> = (q{r, x)T M q{r, x)(q(0, 0)T M q{0, 0))+) 
= (q a (r, x)q /3 (T, x)q x (0, 0)q fi {0, 0)) x 

x(r A/ ) Q/3 ( 7 °r] Vf 7 )A M (6) 

diagrammatically displayed, after subtracting the disconnected contribution, 
in Fig. (1). 

Since the field operators inside the average over the ensemble are already 
time-ordered (being r > 0) we apply Wick's theorem. From the customary 
definitions (see Appendix A): 

(T T [q p (T,x)q x (0,0)]) = S$ x (T,x) (7) 
(T T [q a (r, x) q ,(0, 0)]> = -S£(-r, -x) 

= S^(P-r,-x), (8) 

S F being the quark propagator, it follows that: 

( J m (t, x)JI(0, 0)) = Tt(T m S f (t, x^TlrfSFiP - r, -a;)) . (9) 
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Hence in Fourier space the mesonic 2 point correlation function reads [4]: 



XM (m,p) = -2N c - ]T J 7 ^Tr[r M 5 F (^ n ,fc)7 rj w7 5 F (^„-^,fc-p)], 



VJU (27T) 



3 r d 3 k 

(10) 

where the mesonic frequency is uoi = 2lnT, while in the quark propagators 
u n = (2n+l)7rT, I and n assuming any integer value. The overall factor 2N C 
(N c being the number of colours) comes from the trace over the light flavours 
and colours. 

If one considers the flavour triplet channel the overall flavour factor Nj = 2 
has to be replaced by 



Tr 



T T 



= 25 ab , (11) 

where r a and r b are Pauli matrices acting in the flavour space. The same pre- 
scription allows one to translate the results obtained in the following for the 
interacting case to the non-singlet channel. 



We can now write down explicitly the free meson propagator in the different 
channels getting: 



1 +2? r d 3 k 



\ [i^up) = -2N C - J J^ Tl l S F( iuJ n, k)S F (iu n - iui, k - p)} (12) 
1 +2? r d 3 k 



P) = 2N C - ^ J ^Tr[ 7 5 ^(^ n ,fc) 7 5 5 F (^ n -^,fc-p)] (13) 
1 + °° r d 3 k 

X v (^,p) = -2iV c - /--^Tr[7^(^ n ,fc)7^^ n -^,fc-p)](14) 

X ^( i u; l ,p) = 2N c - J j^Al^SF^k^.S^-iuuk-p)} , 

(15) 

where in the vector (see also Ref. [4]) and pseudovector channels the Dirac 
matrices have been contracted with the Minkowskian metric tensor g^ u . Note 
that in the massless, non interacting case, since {,Sf,7 5 } = 0, x s — X ps an d 

x v =x pv - 

Next we evaluate the imaginary part of the above quantities, thus obtaining 
the free meson spectral function. We sketch the calculation in the vector case 
and just quote the results for the other channels. To simplify the problem we 
confine ourselves to the p = limit. Yet, even in this limit, we can grasp 
important aspects of the plasma physics once the interaction of the quarks 
with the thermal bath is switched on. Indeed <j(uj, 0) (the zero momentum 
projection of the spectral function) is the quantity which is studied on the 
lattice in order to determine whether a meson can survive as a bound state also 
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in the deconfined phase of QCD [12]. Furthermore, as previously mentioned, 
a v (uj,0) is directly related to the back-to-back dilepton production rate [1]. 
Introducing the four-vectors A = (icu n , k) and B = (iu n — iui, k) one can easily 
perform the trace in Eq. (14) obtaining: 

Tr[ 7 ^(4 + m) 7M (^ + m)\ = -8(A ■ B) + 16m 2 , (16) 

m being the bare quark mass. Thus what is left out to compute is the following 
sum over the Matsubara frequencies 

1 v 8K 2 t + £f)-8^ M + 8m 2 

which is easily done with the help of the formulas quoted in Appendix D. 
Then, after performing a proper analytical continuation letting %lo\ — > u + 
(u) being real), one takes the imaginary part of the obtained result. Note that 
the term in round brackets in the numerator of Eq. (17) leads to a purely 
real result (which, incidentally, would require to be renormalized) and hence 
does not contribute. We thus recover the well-known result for the free vector 
spectral function (for u > 0): 



<r>, 0) = -2N C / (03(1 - n{u/2))5{u - 2E k ) ■ 2 • (j| 



UJ 



§g(, - 2mJl- (^)^tanh(./4T) (-2 - (^) 



(18) 



In the above one recognizes the typical threshold behaviour, the square root 
stemming from the integration over the phase space. 

Globally the results for the different channels can be summarized in the com- 
pact expression: 



a^ e (oj,0) = ^(c-2m)^l- ^) 2 c 2 tanh(c/4T) (a + fc^Q , 

(19) 

where the coefficient (a, b) turn out to be given respectively by (1,-1), (1,0), 
(-2,-1), (-2,3) in the scalar, pseudoscalar, vector and pseudovector channel 2 . 



2 These coefficients do not coincide with the ones of Refs. [3,4], due to our choice 
of working with ordinary Dirac matrices satisfying the anti-commutation relation 
{l^,l v } = 2g MJ/ , g^ v being the Minkowskian metric tensor. 
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Fig. 2. The diagram contributing to the quark self-energy in the HTL approximation. 
The wavy line denotes a gluon. The external quark lines are soft. 
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Fig. 3. Dispersion relations to = iv±(k) corresponding to the quasiparticle poles of 
the HTL fermion propagator. Dimensionless units k/m q and tu/m q are used, m q 
being the thermal quark gap mass (see text). The normal (continuous line) and the 
plasmino (dashed line) modes are displayed together with the light-cone (dotted 
line). 

4 Spectral functions in HTL approximation 

In this section we evaluate the meson spectral function in HTL approximation 
for the pseudoscalar and vector channel. In both cases this amounts to dressing 
the quark propagator with the self-energy of Fig. 2, keeping only the term 
proportional to g 2 T 2 , as discussed in detail, for example, in [2]. In performing 
this task it is useful to express the fermionic propagator thus obtained through 
its spectral function (see (B.5)). The latter, in HTL, turns out to read [1,2]: 



p±(u,k) 




[6(u - u)±) + 5{uo + lu t )] + /3±(w, k)6{k 2 




(20) 
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<0/m q 



Fig. 4. Dimensionless (continuum) spectral function m q ■ (3±(uj,k) for space-like 
momenta at k = m q as a function of io/m q . The maximum of /?_ stems from the 
second zero of the function Re(*Al 1 ), occurring in the space-like region. Owing to 
the large imaginary part of the self-energy, this maximum does not show up as a 
well pronounced peak corresponding to a quasi-particle excitation. 



with 



ml ±uj - k 



I3±(u,k) 



2 



+ 



2'"q k 



2 ' 



2 [k(-u±k)+m*(±l-^]n%&)' 

(21) 

where the two dispersion relations u)±(k) are shown in Fig. (3), while the 
behaviour of /3±(u;,k) is plotted in Fig. (4). In the above m q = gT/y/6 is 
the so-called thermal gap mass of the quark, clearly appearing in Fig. (3) at 
k = 0. The HTL spectral functions in Eq. (20) consist of two pieces: a pole 
term, arising from the zeros of the real part of the denominators of *A± in 
Eq. (B.l), and a continuum term, corresponding to the Landau damping of a 
quark propagating in the bath. In the time-like domain u > k the dispersion 
relation u + (k) solves the equation *A~ 1 (a; + (&;), k) = and corresponds to the 
propagation of a quasi-particle with chirality and helicity eigenvalues of the 
same sign. On the other hand u-(k) satisfies the equation *Al 1 (u-(k), k) = 
and describes the propagation of an excitation, referred to as plasmino, with 
negative helicity over chirality ratio. Note that both these excitations are 
undamped at this level of approximation, since the logarithm in *A± doesn't 
develop any imaginary part in the time-like domain. On the contrary, for 
space-like momenta u < k, the quark spectral function receives contribution 
from the imaginary part stemming from the logarithm contained in Eq. (B.2): 
this gives rise to the term f3± in Eq. (20). 

We now address the meson spectral functions in two specific cases. 
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Fig. 5. Diagram giving the pseudo-scalar meson 2-point function in the HTL approx- 
imation. The interaction vertices are 7 5 matrices, while the dressed lines are HTL 
resummed propagators. No vertex correction is required within this scheme. 

4-1 Pseudoscalar channel 



The pseudoscalar vertex receives no HTL correction as shown in [20]. Hence 
in the pseudoscalar channel, following Ref. [4], one has simply to replace, in 
Eq.(13), the free thermal quark propagators with HTL resummed ones (which 
we shall denote with a star), obtaining: 



1 +00 . j3l 

X * s (iu hP ) = 2N c - £ J Tr[7 5 ^(^ n ,fc) 7 5 ^ n -^,fc-p)]. 



"PJU (27T) 



(22) 



Setting q = k — p and making use of the spectral representation (B.6) of the 
quark propagator one gets: 



1 +£? r d 3 k r+°° r+°° 1 1 
X ps (tuj h p)=2N c - 77^ / du! / du 2 -. , ; X 

Pn^Z^J {27r) 6 J -00 J -00 lU) n — U)\ lLU n — lUJl — UJ 2 

x Tr[ 7 5 pHT L (^i, fc)7 5 PHTL(^2, q)\ ■ (23) 

Summing over the Matsubara frequencies with a standard contour integration, 
performing the usual analytical continuation iu\ — > uj + ir] + and taking the 
imaginary part of the result thus obtained one gets: 



aP s (cu,p) 



-2N r 



d?k r+00 r+00 

——-(e^-l) du x duj 2 h(uJi)n(uJ2)~x 

[Z7l)° J-00 J-00 



x 5(u -uj!- lo 2 ) ■ Tr[7 5 p H TL(^i, fc)7 5 PHTL(-^2, q)\ ■ (24) 
Now, inserting Eq. (B.5) into Eq. (24) and since 



Tr 
Tr 



7 ° Tl . k 5 7 T-7 • q 
7 n 7 o 



7O Tl . k 5 7 ° ±7 • q 
7 n 7 „ 



-(1-k-q) 
-(1 + k-q) 



(25a) 
(25b) 
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one gets [4]: 



f U K r+oo r+oo 

a ps (uj,p)=2N c J J duj x J <Lj 2 n{uj 1 )n{uj 2 )8{uj-u x ^j 2 )x 

x {(1 + fc • g)[/o+(wi,A;)/) + (w2,g) + /o_(wi,A;)/)_(w2,g)]+ 

+ (1 - k • q)[p + (uJi, k)p„(uj 2 ,q) +p-(wi, k)p + (u 2 ,q)]j . (26) 

where the identity p + (—u, k) = p-.(u, k) has been used. 
In the case p = the above formula reduces to: 



<7 pB (a;,0) = — -(e^-1) / rffcA; 2 / dui du 2 n(ui)n(ui) 



7T 2 ' JO J-oc 



5(o; - tui - ^ 2 )[p+(^i, k)p + (uj 2 , k) + /0_(o>i, k)p-(oj 2 , k)] . 



(27) 



Inserting then into Eq. (27) the explicit expression for p± given in Eq. (20), 
one finds, as first pointed out in Ref. [1], that the meson spectral function 
consists of the sum of three terms: pole-pole (pp), pole-cut (pc) and cut-cut 
fee). These different contributions will be discussed in Sec. 6. 



4-2 Vector channel 




Fig. 6. The diagrams contributing to the vector meson spectral function in the HTL 
approximation. Dressed lines and vertices denote HTL resummed propagators and 
HTL effective quark-photon vertices respectively. The momentum flowing in the 
loops dressing the vertices is hard. The external dashed lines denote a photon (or a 
probe with the same quantum- numbers). 

The vector meson spectral function in the deconfined phase, being related to 
the thermal dilepton production rate from the plasma, is clearly of importance 
for the QGP phenomenology. When computed in the HTL approximation it 
should provide a reliable estimate for the soft dilepton production rate. For 
the sake of completeness we quote here the HTL vector spectral function 
previously obtained in [1,2,3,4], whose meaning in term of Feynman diagrams 
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is illustrated in Fig. (6). Its expression is considerably more cumbersome than 
the pseudoscalar one reading [1,2,3,4] 



<7 v (a;,0) = -^V" 



x 



4 1 



IT' 



r+oo 

■1) / c 
Jo 



dkk 2 



+oo 



+oo 



dui I duj 2 n{uji)n{uj 2 )8{uj—uji—uj 2 ) 



2ku 



p+(ui, k)p-(u 2 , k) 



1 + 



uo( + u£ - 2k 2 - 2m\ 



2\ 2 



2ku 



x p+(ui, k)p + (u 2 , k) + 1 - 



- 2k 2 - 2m\ 
2ku 



2\ 2 



p_(^i, k)p-(u 2 , k) 



+Q{k 2 -ujl) 



4Axj 2 



(l + ^WK*) + (l-y)P-K*) 



(28) 



and is derived by dressing in Eq. (14) both the propagation lines (replacing 
the free thermal quark propagators with HTL resummed ones) and the bare 
vertices. For the latter the following effective HTL vertex for the coupling of 
a quark with a vector probe is adopted: 



r iI (P l ,P 2 )= lfl + m 2 q J 



4:71 



(P 1 -K)(P 2 -K) 



(29) 



where K — (l,k), P\ = (icu ni , pi) is the four-momentum of the incoming quark, 
while P 2 = (iu n , 2 ,p 2 ) refers to the outgoing one. Again we note that, because 
our 7^ matrices satisfy the usual Dirac algebra (with the Minkowskian metric 
tensor g^ u ), we get a different sign in front of the vertex correction with re- 
spect to what quoted in Ref. [2]. 



At this point a careful analysis of the momenta over which one performs 
the integration is appropriate. In fact, according to the original proposal by 
Braaten and Pisarski [21] for constructing an effective theory for the soft modes 
in the QGP, one has to replace every soft fermionic line (i.e. with K ~ gT) with 
an HTL resummed propagator, while the hard propagation lines should be left 
undressed. This is because, for K ~ gT, also the self-energy Ehtl(-^) ~ gT: 
hence all the diagrams with repeated self-energy insertions in the fermionic 
lines give contributions of the same order, i.e. 

zi . (-s)Zi . (-£)— ^ — . 

If If. K If gT 

In addition, if all the lines entering into a vertex are soft one has to replace 
it with the HTL one given by Eq. (29), otherwise the vertex should be left 
untouched. On the other hand a hard fermion propagator (i.e. with K ~ T) 
can be expanded perturbatively, since the kinetic term dominates over the 
interaction, and in first approximation one can keep its non-interacting ex- 
pression. Hence the integration over k in Eq. (28) really should cover only 
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the soft regime of momenta, while in the integration over the hard momenta 
domain one should use the bare quantities. 

Actually the vector meson spectral function in the deconfmed phase and the 
associated thermal dilepton rate have been recently evaluated on the lattice 
[22,23,24]. It turns out that their behaviour at low energies is very different 
from the HTL result given by Eq. (28). Specifically, whereas the HTL spec- 
tral function develops an infrared divergence for oj — > (stemming from the 
cut-cut contribution) [1], the lattice result drops to zero below uo/T ~ 3, sug- 
gesting a threshold effect due to a thermal mass acquired by the quarks. In 
this connection we remind the reader that the spectral function of a meson is 
not measured directly on the lattice, but has to be obtained by inverting Eq. 
(5). A procedure referred to as Maximum Entropy Method [12,13,14,22,23,24] 
is usually employed to accomplish this task. 

Finally, for what concern the experimental results, in a heavy ion collision 
dilepton pairs are produced also in processes different from the annihilation 
of qq quasi-particles in the plasma phase 3 . To get an estimate for the produc- 
tion rate which can be eventually compared with the experiment, one has to 
account also for other sources [5] : in particular the (non-thermal) contribution 
due to Drell-Yan processes, the dileptons coming from the hadronic phase and 
the decay of vector mesons after the freeze-out. 

In the following we shall concentrate on the pseudoscalar channel only, aiming 
at an improved treatment of the latter. 



5 Beyond HTL 

In this section we attempt to improve upon the quark spectral functions ap- 
pearing in the loop of Fig. (5). For this purpose we observe again that, in 
making the convolution of the two fermionic propagators, one has to integrate 
over all the scale of momenta (hard and soft). Now, while the HTL approxi- 
mation is supposed to dress correctly the propagation of the soft modes, this 
is not so for what concerns the hard modes. Thus, by replacing naively the free 
thermal quark propagators in Eq. (10) with HTL resummed ones, one treats 

3 Actually, even in the deconfmed phase as we will discuss in detail for the pseu- 
doscalar case, the HTL meson spectral function gets contributions from far from 
trivial many-body processes. 
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(a) (b) 

Fig. 7. Next to Leading corrections to a hard quark propagator, arising from the 
interaction with a soft transverse (a) and longitudinal (b) gluon: for these HTL 
resummed propagators are used. Thus the two diagrams receive contributions from 
all orders in perturbation theory. 

incorrectly the contribution to the integral arising from hard momenta. This 
is analogous to what happens when one tries to evaluate thermodynamical 
quantities for the QGP, like the entropy and the baryon density, in a pure 
HTL approximation. As pointed out in Ref. [7], in HTL one gets the right 
contribution of order g 2 to such quantities, but only part of the g 3 term (actu- 
ally this one is strictly a non-perturbative contribution, being non analytical 
in a s = g 2 /47r), namely the contribution arising from the soft modes. 



Now we shortly examine how the Next to Leading contributions, amounting to 
dress properly the propagation of the hard degrees of freedom in a QGP, are 
kept into account in a self-consistent calculation of quantities like the entropy, 
the baryon density and the quark-number susceptibilities [7,8,10]. 
These observables receive a non-analytical contribution of order g 3 both from 
the soft modes and from the proper dressing of the hard propagators. To get 
the latter one has to evaluate the correction 5Ti±(u = ±fc) 4 to the self energy 
of a normal quark mode (with hard momentum). Such a contribution corre- 
sponds to the two diagrams displayed in Fig. 7 where the hard quark interacts 
with a soft gluon, longitudinal and transverse, respectively, whose resummed 
propagators are the HTL ones (being the gluons soft such an approximation 
is fully justified). Note that this is a non-perturbative procedure, the two di- 
agrams in Fig. 7 getting contributions from all orders of perturbation theory. 
Actually the correction 8H±{u = ±k) is a far from trivial function of the 
momentum and generally it can be dealt with only numerically. Nevertheless 
its contribution of order g 3 to the entropy and the baryon density can be 
analytically evaluated. Indeed in [7,8,10] the authors replace this momentum- 
dependent self-energy correction with a constant shift of the asymptotic quark 

4 Note that only the light-cone projection of the hard quark self-energy enters into 
the thermodynamical calculations. The same is true for the hard transverse gluon 
as well. Since on the light-cone the HTL self-energies coincide with the full one-loop 
calculation, the entropy and the baryon density might be privileged quantities to 
study within the HTL approximation. 
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mass in order to reproduce the right coefficients of the terms of order g 3 com- 
ing from the hard sector for the entropy and the baryon density 5 . 
This mass correction turns out to be negative, reading [7,8,10]: 

Sm 2 00 = -^-g 2 C f Tm D , (30) 

where (in the case of zero chemical potential) the Debye screening mass in the 
HTL approximation is: 



m D 



being Nf the number of flavours and C/ = 4/3 the quadratic Casimir operator 
for the fundamental representation of 577(3). 

It is worth stressing that at the experimentally accessible temperatures (i.e. a 
few times T c ) g > 1, hence the (negative) correction (30) becomes of the same 
order of magnitude of the HTL asymptotic mass which reads 

ml = 2m 2 q = ^ , (32) 
thus possibly leading to a tachyonic singularity. 

This problem was addressed in Ref. [7] where two possible solutions where 
proposed. The first amounts to replace the expression ml = ml + 5ml with 
its first Pade' approximant, namely 

ml = ml + 6ml — < = — %r ■ ( 33 ) 

m l 



The other starts with the correct equation for the Next to Leading quark 
asymptotic mass m^, namely 



ml=m 2 QO -—g 2 C f ^N c + ^-Tm 00 , (34) 

and replaces — > in the term linear in m^, thus obtaining the "self- 
consistent like" equation: 



ml = ml - — g 2 C f ^ N c + -f T , (35) 

whose positive root provides the asymptotic mass for the hard quark. In the 
following, in accord with the literature [7], we shall refer to such an approach 

5 Analogously the correction SW.t (lo = ±k) to the self-energy of a hard momentum 
transverse gluon is absorbed in a shift of its asymptotic mass. 
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as the Next to Leading Approximation (NLA). 

To replace in the term linear in rh^ the HTL mass with its "self-consistent" 
value may appear somewhat arbitrary. Yet the effects of such a procedure were 
tested in a scalar toy model [7] , where one can solve the gap equation exactly 
(although numerically), and the discrepancy between the exact result and the 
above discussed procedure turned out to be ~ 1% even for large values of the 
coupling (g ~ 2), at least as far as the entropy is concerned. 
The NLA asymptotic mass applies to the hard quark modes. Since the 
soft modes are correctly described by the HTL approximation, the strategy 
followed in Refs. [7,8,10] was to single them out by introducing a cutoff A at 
an intermediate scale of momenta. A good choice for the latter is: 



which represents the geometric mean of the spacing between the Matsubara 
frequencies (hard scale) and the Debye screening mass in HTL (soft scale). 
Later we will investigate the sensitivity of the results to the value of this cut- 



Now, in accord with the above discussion, in the integration over the spatial 
momenta in the formulas for the entropy and the baryon density one is allowed 
to keep the HTL expression of the integrand for k < A (k being the modulus 
of the three-momentum). Hence the entropy will receive contribution from the 
quasi-particle poles (normal quark mode, plasmino, transverse and longitudi- 
nal gluons) and from the Landau damping as well. For k > A instead one 
takes into account only the propagation of the physical modes (normal quarks 
and transverse gluons), which are now dressed by the asymptotic masses ob- 
tained in NLA, neglecting the damping contribution. A comparisons of the 
results obtained in such a framework with the one coming from lattice cal- 
culations are reported in Refs. [7,8,10] and the agreement appears satisfactory. 

A similar strategy is here applied to the evaluation of mesonic correlation 
functions. In particular we explore in detail how an improved approximation 
at the level of the quark propagator is reflected on the (zero momentum) 
meson spectral function (in the pseudo-scalar channel) c(uj,p = 0) and thus 
on the temporal correlation function G(r,p = 0) obtained through Eq. (5). 
Before addressing the issue of improving the HTL results, let us discuss the 
basic features of the HTL quark spectral function that we want to keep in the 
new scheme. 

The structure of the spectral representation given in Eq. (B.3) guarantees 
that: 



*S being the HTL resummed quark propagator. This expresses the fact that, in 
spite of the appearance in the spectrum of a gap mass m q and of an asymptotic 
mass rhoo for the normal quark mode, both temperature dependent, neverthe- 




(36) 



off. 



{^,7 5 } = 



(37) 
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less the chirality of a quark propagating in the thermal bath is conserved. 
One would like Eq. (37) to hold even when a better description for the hard 
fermionic modes is given. Thus, in analogy with Eqs. (B.5) and (B.6), the 
following spectral representation for the NLA quark propagator 



S NLA (tuj n , k) = - [ + °°du pNLA{u; ' k) , (38) 



o ytu n , K) = - I 

where 



P NLA (-, fe) = ^^pf A (-> k) + ±±£±^{u,, k) , (39) 

is assumed to hold. The functions p± LA are discussed below. 

Now the HTL quark spectral functions p± obey the following sum rule [2] 

/+oo 
dwp±(u,k) = 1 , (40) 
-oo 



which can be cast in the form 

r+k 
-k 

where the functions 



r+k 

Z + (k) + Z_(k) + dwp±(ou,k) = 1 , (41) 

J —k 



are the residues of the quasi-particle excitations. 

Since Eq. (40) is derived directly from the equal-time anticommutation rela- 
tion: 

ty a (f , x), ^(t, x')} = 5 a(3 5(x - x') (43) 
one would like it to hold for the NLA quark spectral function as well. 
It is of interest to discuss how the saturation of the sum rule (40) is im- 
plemented in the HTL approximation. To this aim we recall the asymptotic 
behaviour of the HTL fermionic dispersion relations for k ^> m q which is found 
to be 



u+(k) ~ \jk 2 + m 2 ^ , u-(k) ~ k + 2k exp ^- 



2k 2 + m 2 q 



where the HTL quark asymptotic mass has been defined in Eq. (32). Concern- 
ing the residues one gets, always for k ^> m q , 

From the above it follows that: 

limZ + (Jfe) = l and lim Z_(k) = . (46) 

fc^oo k— >oo 



17 



Hence, for k ^> m q (admittedly a situation where the HTL approximation is 
no longer realistic) the normal quark mode tends to saturate the sum rule (41) 
and its dispersion relation approaches the one of a free particle of mass m^. 



In line with the above arguments we propose an improved expression for the 
quark spectral function which satisfies the following constraints: 

• for small momenta it reduces to the HTL one; 

• for large momenta it yields a spectrum dominated by an undamped exci- 
tation whose dispersion relation approaches e^ LA = \Jk 2 + m^, being 
the solution of Eq. (35); 

• it obeys the sum rule (40); 

• the associated quark propagator anti-commutes with 7 5 . 

A prescription should also be given to allow a transition from the soft to the 
hard momenta regime. In this connection, keeping for the intermediate cutoff 
A the value given in Eq. (36), a natural choice might appear the following one: 

pl LA (u, k) = 9(A - k)p±(u>, k) + 0(k - \)6(u =F e£ L A ) , (47) 



where p± are the HTL quark spectral functions, while e^ LA = \Jk 2 + 
being the NLA quark asymptotic mass. 
Actually, since u;+(A) > e^ LA (A; = A), a narrow interval of values of u (i.e. 
2e NLA (A) < uj < 2u + (A)) exists in which the meson spectral function would 
get contribution both from the normal quark mode u + (k), namely 

2N r A 

-^-{e^ - 1) / dk k 2 n 2 {u + {k))Z 2 Jk)5{uj - 2uo + {k)) , (48) 

7T Z JO 

and from the asymptotic quark dispersion relation e^ LA , namely 

^V" - 1) rdkk 2 n\e^)5{u - 2e f A ) • (49) 

Hence, due to this overcounting of the same physical mode, the meson spectral 
function would display an unphysical bump for 2e NLA (A) <u< 2cj + (A), cr(u) 
being discontinuous at the border of the interval. 

In order to cure this problem (at least partially, as we will see in the following) 
and to interpolate smoothly between the soft and hard regimes we introduce 
two additional momenta Ai and A 2 defined as follows: 



w + (A!) = y /A^ + ml (50a) 
c + (A) = v /A| + m^ . (50b) 

Hence the following guess on the spectral functions p^ LA (u;, k) entering into 
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(a) (b) 

Fig. 8. Higher order diagrams contributing to the pseudoscalar meson spectral func- 
tion. All the fermionic lines are hard. The dressed wavy lines are HTL resummed 
gluon propagators. The effect of diagram (a) (together with the corresponding one 
with the self energy insertion in the other fermionic line) are taken into account 
in an effective way through the introduction of the NLA quark effective mass. The 
effects of diagram (b) , which in principle might give a contribution of the same order 
of (a), are not accounted for in our NLA scheme. 



Eq. (39), namely 



pl LA (cu,k) 



with 



p±(u, k) if k < A 1 

cos 2 (a(k))p±(uo, k) + sin 2 (a(k))5(uj T e^ LA ) if Ai < k < A 2 



8{u T e% LA ) 



if k > A 2 



a{k) = | 



fc-Ai 
A 2 -Ai 



(51) 
(52) 



appears suitable. 

Replacing then the HTL quark spectral functions entering into Eq. (27) with 
the above given ansatz, one gets the following NLA expression for the pseudo- 
scalar meson spectral function (for uj > 0): 



<7mt a (uj, 



2N C 



+00 



dk k z 



duj\ 



+00 



duj%n(uj{\n(u}-i, 



NLAl^-/ _ 2 

8(uj-uj x -uj 2 )\p^(uj x , k)pl LA (uj 2 , k) + pT^(ui, k)p^ A (u 2 , k)} . (53) 



7T JO J —00 J —00 

NLA/ . 1 \ NLA/ . i\ , NLA/ . i \ NLA/ 



In conformity with Eq. (5), one can next obtain the associated zero momen- 
tum temporal correlator G^ la (t, 0) whose behaviour will be later shown for 
different values of the temperature. 



Before presenting our numerical results we stress that our approach is just an 
effective way of accounting for the contribution to the meson correlation func- 
tion arising from the hard fermionic modes, for which the HTL approximation 
no longer holds. Specifically we embody, in analogy with entropy calculations 
of Refs. [7,8,9,10], the effects of diagram (a) in Fig. (8) into the NLA quark 
asymptotic mass, which turns out to be lower than the HTL one. On the other 
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hand we neglect the corrections arising from diagram (b) which in principle 
could be equally important. 



6 Numerical results 



In this section we first assess the weight of the different contributions (pole- 
pole, pole-cut and cut-cut) to the HTL pseudoscalar spectral function 

a^ TL (uj, 0) = a pp (uj, 0) + a pc (uj, 0) + a cc (to, 0) . (54) 

Next we explore the modifications of the HTL results stemming from the more 
proper NLA treatment of the hard modes. 

For the sake of consistency, we start from the explicit expression for the three 
terms involved in a^ Th [ui,0). These have already been deduced in Ref. [4]. 
For zero spatial momentum the pole-pole term, which selects the quasi-particle 
peaks in both quark spectral functions in Eq. (27), is given by: 



at e fa> _ i 



2tt 2 m\ 



+ 2EnK(fc 2 ))[l-^fe))](^fe)-^)(^fe)-fc 2 2 ) l „, n. ^ 1 71 m + 

k 2 



iv' + (k 2 )-u;'_(k 2 )\ 
k 2 



k 3 Z\uj_(k 3 )\ 



(55) 



In the above k\ is the unique solution of the equation u — 2u + (ki) = 0, whereas 
two are the values of k 2 fulfilling u — uj + {k 2 ) + uj^{k 2 ) = and, likewise, two 
values of k$ satisfy oo — 2uo-{k 3 ) = 0. 

In order to identify the different contributions to the pole-pole term with well 
defined physical processes let us introduce the notation q + (q + ) for the normal 
(anti-) quark mode, g_ (g_) for the (anti-) plasmino mode and M for the 
pseudo-scalar meson. 

Thus the first term in Eq. (55) corresponds to the annihilation of a normal 
quark anti-quark pair into a meson at rest 

Q+ + Q+ — > M , 

which clearly turns out to be proportional to the mean occupation number of 
normal (anti-) quarks states, hence the factor h 2 (uj + (k 1 )). Due to energy con- 
servation the density of states turns out to be proportional to \d(2uj + ) / dk]' 1 . 
On the other hand the second term is related to the decay of a normal quark 
mode into a plasmino mode with the same momentum and a meson at rest, 
namely 

q+ >q- + M . 
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Therefore in this term one of the statistical factor accounts for the mean 
occupation number of normal quarks states (hence the Fermi distribution 
n(u + (k2))), while the other one expresses the fact that the final plasmino 
has to occupy empty states (hence the Pauli-blocking factor 1 — n(u-(k2)))- 
In this case the density of states is proportional to \d(uo + — uo-)/dk\~ l . 
Finally the third term describes the annihilation of a plasmino anti-plasmino 
pair into a meson 

q- + q- > M , 

whence the statistical n 2 (o>_(/c 3 )) and the density of states \d(2uo^)/dk\^ 1 fac- 
tors. 

It is important to realize that, owing to the minimum in the plasmino dis- 
persion relation, two values of uo exist for which the denominators in Eq. (55) 
^+(^2) — (^2)1 and |u/_(/j3)| vanish, thus yielding singularities in the meson 
spectral function . These are the so-called Van Hove singularities, essentially 
reflecting a divergence in the density of states. One finds (numerically) that 
they occur at uo m 0A7m q and at uo ss 1.86m g . The impact of these singu- 
larities on the meson spectral function is dramatic for soft values of uo. Here 
indeed the HTL meson spectral function turns out to be orders of magnitude 
larger than the free one, as can be seen in Fig. (9). 

Turning to the pole-cut contribution (for positive energies) it is found to be: 



a pc (uo) 



2N r - 1 



71 z 



mi 



f 



dk k 



6{k 2 — {uo — uo + ) 2 )n(uo — uo + )n{uo + )l3 + {uo — uo + , k)(u+ — k 2 ) 
+6(k 2 -(uo- uo^) 2 )n(uo - oo_)n{oo-)(3-{oo - uo^, k)(uo 2 _ - k 2 ) 



(56) 



Finally the cut-cut contribution is given by: 



a cc (uo) 



2N r 



7T^ 



-1) 



dk k z 



dx n(x)n(uo — x)6(k 2 — {uo — x) 2 )- 



\P+(x, k)j3 + (uo — x, k) + P-(x, k)j3_(uo — x, 



(57) 



The three above quoted contributions to the meson spectral function in the 
pseudoscalar channel in a pure HTL approximation are plotted in Fig. 9: they 
numerically coincide with the results previously obtained by Karsch et al. [4]. 



The question we now raise is whether the behaviour of the meson spectral 
function at soft energies, displayed in Fig. 9, is qualitatively modified once a 
better description of the hard fermionic modes is adopted. 
Since the divergence in the density of states, which gives rise to the Van Hove 
singularities in the pole-pole term, stems from the stationary points of the 
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Fig. 9. The various contributions (pole-pole, pole-cut and cut-cut) to the dimension- 
less spectral function of a pseudoscalar meson a ps /T 2 at zero spatial momentum as 
a function of x = u/T. Here the quark propagators have been evaluated in the HTL 
approximation both for hard and soft momenta. In order to test the present results 
versus the ones obtained in Ref. [4], here the plot is given for a value of T such 
that g(T) = entailing m q = T. In the pole-pole contribution one recognizes the 
occurrence of the Van Hove singularities at x = 0.47 and at x = 1.856. Also plotted 
is the free spectral function. 



curves given by the equations 



Lu = 2tu_(k) (58a) 
and 

u = iu + (k) -uj_{k) , (58b) 

to answer the above question one should check whether the choice of the inter- 
mediate cutoff A separating the soft modes from the hard ones yields A > ki, 
where ki are the stationary points of the curves (58a) and (58b). In such a 
case indeed the Van Hove singularities will still be present in the spectral func- 
tion. Therefore we report below the numerical values of A 1; A, A 2 (normalized 
to the quark thermal mass m q ) for two different choices of the intermediate 
cutoff A, for T = 2T C and T = 4T C . We also report the corresponding value of 
the gauge coupling g obtained through the QCD /3-function as discussed below. 



T = 2T C (0=1.785) 


A 1 /m q 


A/m q 


A 2 /m q 


A = y/nTrha 


3.312 


3.492 


3.665 


A = ^2nTm d 


4.806 


4.939 


5.069 
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T = AT C (0=1.548) 


A 1 /m q 


A/m q 


A 2 /m q 


A = y/nTrhd 


3.593 


3.750 


3.902 


A = V2vrTm d 


5.188 


5.303 


5.417 



From the comparison with the curves in Figs. 3 and 10 which, being normalized 
to the quark thermal mass m q , hold for every value of the temperature, one 
realizes that the stationary points leading to the Van Hove singularities always 
occur in the soft domain k < Ai. Hence we conclude that, qualitatively, the 
behaviour of the HTL meson spectral function does not change once a better 
treatment of the hard fermionic modes is kept into account. The quantitative 
changes are below examined. 




k/m 



Fig. 10. The behaviour of the difference between the two fermionic dispersion re- 
lations Cj + (x) — uj^(x) . The dimensionless variables uj = uj/m q and x = k/m q are 
used. The maximum of the curve occurs at x = 1.179 and corresponds to a value of 
Co = 0.4698. 

Accordingly, we next present the results obtained in the Next to Leading Ap- 
proximation (NLA), which has been briefly discussed in the previous section 
and which allows a better treatment of the contribution of the hard quark 
modes in the meson spectral function. 

Before addressing the discussion of the numerical results a few words, regard- 
ing how the running of a s is related to the temperature, are in order. 
The running of the QCD coupling g with the renormalization scale p, (intro- 
duced in the dimensional regularization scheme) is here obtained through the 
two- loop /9-function with Nf = 2, in accord with [7]. The modified Minimal 
Subtraction (MS) renormalization scheme is employed. Dealing with massless 
quarks, the only physical scale entering into the problem is the temperature. 
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Hence the evaluation of the coupling g at the scale ft = 2ttT, which corre- 
sponds to the spacing between the Matsubara frequencies, appears a reason- 
able choice. For further details see [7] and references therein. 



In Fig. 11 we display the ratio of the Pade' and NLA asymptotic mass with 
respect to its HTL value. The two curves represent the solution of Eqs. (33) 
and (35). It appears clearly that the interaction of the hard quark modes with 
soft longitudinal and transverse gluons (for which one uses HTL resummed 
propagators) lowers the value of the asymptotic quark mass and that the two 
approximation schemes lead to similar results. 

In Fig. 12 the behaviour of the ratio A/T as a function of the temperature is 
given. 



: NLA 

: |-- ■ Pade' 

0.8 - 



0.6 - 




0.2 - 



C , 1 , 1 , 1 , 1 , 1 , 3 

2 2.5 3 3.5 4 4.5 5 

T/T 

c 

Fig. 11. Ratio of the asymptotic quark mass over its HTL value as a function of 
T/T c : here the (negative) correction Sm^ has been included in two different ways: 
with the first Pade' approximant (dot-dashed line) and in the "self-consistent like" 
way (continuous line). 

In Fig. 13 we plot the pseudoscalar meson spectral function obtained in the 
NLA scheme with two different choices of the cut-off parameter A. This is 
compared with the pure HTL result, with the curve obtained assigning to the 
quarks a mass m^ jA and with the non interacting case. Notice that the curves 
NLAl and NLA2 display a peak which has no physical meaning, but is an 
artifact of our approximation for the quark spectral function. In fact a narrow 
energy interval exists where the meson spectral function gets contribution both 
from the normal HTL quark mode and from the NLA asymptotic dispersion 
relation. On the other hand, as discussed in the previous section, without the 
prescription (51) the spectral function would display an even worse unphysi- 
cal behaviour. Although our findings for this intermediate scale of frequencies 
cannot be trusted, nevertheless for the soft and hard regimes they are reliable. 



24 




T/T 

c 



Fig. 12. Temperature behaviour of the intermediate cutoff A = \/2TxTmD separat- 
ing soft and hard momenta. Dimensionless units A/T and T/T c are adopted. The 
decreasing of the ratio A/T reflects the running of the strong coupling constant. 

In Fig. 13 the kernel entering into the definition of the meson temporal corre- 
lation function (see formula (5)), namely 

K(t, w) = cosh(cu(r - (3/2))/ nnh(wP/2) , (59) 

is also plotted for r = (5/2. One infers from the behaviour of K{r = (3/2, cu) 
that, for this value of r, the contribution from the soft sector to the integral 
in Eq. (5) is significant. 

Finally in Figs. 14 and 15 we show, for two different values of the temperature, 
the zero-momentum temporal correlator G(t = —ir,p = 0) in the two differ- 
ent approximations schemes discussed in this paper: HTL and NLA. We also 
show the curve obtained dressing the quark modes with the NLA asymptotic 
mass and setting to zero the plasmino and the damping contribution to the 
quark spectral function. 

As expected from the general theory, the curves are symmetric for r — > (3 — r. 
The correlation function diverges for r — > {(3) independently from the ap- 
proximation scheme employed and even in the non-interacting case. Hence 
the meaningful quantity to study in order to assess the impact of the inter- 
action with the thermal bath on the meson correlation function is the ratio 
G , (r)/G frcc (r) given in panel B in Figs. 14 and 15. Note that for r — >• 0(/3) one 
recovers the free result. 

For r sufficiently small (or close to (3) the NLA curves stay above the HTL 
result. This may be understood on the following basis. As r gets smaller, the 
contribution to G(t) coming from large frequencies grows, as one can see from 
Eq. (5). But for large values of u> the meson spectral function (both in HTL and 
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Fig. 13. Behaviour of the zero momentum pseudoscalar spectral function a ps /T^ 
versus u>/T c in different approximations: free result, HTL, NLA and quarks endowed 
with a thermal mass ttinla- NLA1 corresponds to the choice A = v / 2vrTrfiD, NLA2 
to A = ^irTm D . The plot refers to T = 2T C . 




Fig. 14. (a): Behaviour of G(r)/T c 3 vs t//3. (b): Behaviour of G(r)/G free (r) vs t//3. 
NLA1 corresponds to A = \/2TrTmo, NLA2 corresponds to A = \/irTifiE>. In panel 
(b) we also display the result obtained in the case of quarks endowed with a thermal 
mass m = ttinla- The curves are given for T = 2T C . 
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Fig. 15. The same as in Fig. 14, but for T = 4T C . 

in NLA) is dominated by the normal quark modes, in which mj^ LA ) < m^ TL ) . 
On the contrary, for r close to (3/2, where the contribution from soft frequen- 
cies becomes more important, the HTL result stays significantly above the 
NLA curves. 

It is also evident that the NLA curve are not much sensitive to the choice of 
the cutoff A and that, as one moves to larger values of the temperature, one 
approaches (slowly) the free result. 

In any case, in the whole range < r < (3, is G(t) < G frcc (r). This outcome 
appears to agree with the results for the correlation function along the z-axis 
obtained in [25]. Here, on the basis of a dimensionally reduced effective la- 
grangian, a positive correction (of the order of 5%) was found to occur for the 
meson screening mass. 



7 Conclusions 



In this paper we have explored the predictions of the NLA framework on the 
thermal meson correlation functions in the deconfmed phase of QCD, thus 
extending past work of Blaizot et al. who successfully investigated, employ- 
ing the NLA scheme, the thermodynamical properties of the QGP. Notably 
their results turned out to be well supported by lattice calculations, even in 
a regime of temperatures where g is not small (g ~ 1 or larger), in which 
the separation of different momentum scales (i.e. g 2 T, gT, T) might appear 
questionable. Hence we felt justified to employ here the same scheme to deal 
with other observables, like the meson correlation functions. 
Actually what we have proposed is a variant of the NLA scheme, smoothly 
matching the soft and hard momenta regimes. 

With a more careful treatment of the contribution to the meson correlation 
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function arising from the hard quark modes we have still found that the NLA 
does not change dramatically the HTL results which turns out to be quite 
robust. In particular the peculiar behaviour of the meson spectral function for 
soft energies (Van-Hove singularities) remains (almost) unaltered. 
The detection of the V.H. singularities in the meson spectral functions, be- 
ing related to the minimum in the plasmino dispersion relation, would be a 
clear signal of deconfmement. Indeed the doubling of the fermionic disper- 
sion relation for non-vanishing three-momentum is not strictly related to the 
HTL approximation, whose validity, for values of the temperature just above 
T c , is questionable. Actually, as discussed in [26,27], it is possible to show, in 
a totally non perturbative framework, that in the high temperature chirally 
symmetric phase of QCD the fermion propagator does always display two dis- 
tinct dispersion relations, one of which developing a minimum. 
Unfortunately the possibility of getting an experimental evidence of these sin- 
gularities, for instance from the dilepton production rate for what concerns 
the vector channel, is dim, owing to the fact that during the expansion of 
the fireball their position moves toward lower energies as the temperature de- 
creases (experiments provide time-integrated yields). Furthermore there are 
other sources of dileptons which may cover the signal. 

At present, lattice calculations of (zero momentum) meson spectral functions 
(of light quarks) <7m(o;) do not show any evidence of sharp resonances in the 
soft energy domain which could be related to the V.H. singularities [22,23,24]. 
This occurrence may depend on the fact that what one actually measures on 
the lattice is the temporal correlator Gm(j) on a finite set of euclidean times. 
Hence, to get ct m (uj) from the lattice requires to invert Eq. (5) and the left 
hand side of this equation is known only on a finite set of values of r. To per- 
form this inversion, the Maximum Entropy Method has been recently adopted 
[12,13,14,22,23,24], which in any case requires some a priori information both 
on Gm(j) and on a M {uj). As a consequence it is not possible at present to as- 
sure that the spectral function reconstructed in this way is unique and that the 
existence of very sharp resonances like the V.H. singularities can be brought 
to evidence with this method. 

Also of interest is to point out that the V.H. singularities, associated to the 
vanishing of the group velocity of a propagating wave, also occur in other con- 
tests in many-body physics. For example in a non relativistic ferromagnetic 
system at T = undergoing a spontaneous symmetry breaking {Quantum 
Phase Transition), the dispersion relation of the Goldstone mode (magnon) 
displays a maximum, at which the group velocity is zero [28]. 
A promising route to the study of mesonic correlations in the deconfined phase 
of QCD is offered by the possibility of applying the same strategy developed 
in this paper to the evaluation of meson screening masses (i.e. considering the 
propagation along a spatial direction) and comparing the results with those 
obtained in different approaches (for instance on the lattice or with effective 
lagrangians) . Indeed, due to the larger number of sites along the spatial direc- 
tions with respect to the (euclidean) temporal one, lattice calculations should 
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provide more reliable results for this kind of observables. 
Concerning the shortcoming of the framework employed in this paper, it should 
be pointed out that we took into account corrections to the free meson corre- 
lation functions arising merely from dressing the fermionic propagators, due 
to their interactions with the medium (i.e. the gluons of the thermal bath). 
On the other hand, if, within the HTL scheme, one is allowed to leave the 
pseudoscalar vertex untouched (since the HTL correction to the pseudoscalar 
vertex vanishes), this is no longer true when one tries to go beyond the HTL 
approximation. While we have not accounted for this fact in our analysis, it 
certainly represents an item worth to be addressed. 

Equally important appears to be the study of the chiral and of the quark 
number susceptibilities in a scheme improving upon the HTL results obtained 
in Refs. [29,30]. 

Finally we stress that we don't expect the framework employed in this pa- 
per to give reliable results for temperature values lower than 2.5 -i- 3T C . For 
the range of T just above T c other schemes [15,16], predicting the presence of 
bound states also in the deconfmed phase, appear more realistic. 
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A Spectral representation of the free quark propagator 

In configuration space the Matsubara fermionic propagator is defined as fol- 
lows: 



whereas in Fourier space, for the non interacting Dirac field, it reads: 



S afS (T, x; t', x') = (T T (ip a (-ir, x)i\) & {-%t' , x'))) 

= 9(T-r'){ij a (-iT,x)M-iT',x')) 
-0(t' - r){ip p {-iT J ,x')i/> a (-iT,x)) , 



(A.l) 



S F (iuj n ,p) 



$ + m 



p2 _ m 2 

jf + m 



(A.2) 
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In the above po = iw n , E 2 = p 2 + m 2 and, as in Blaizot and Iancu in [11], the 
matrices 7 M (fi = 0, . . . 3) satisfy the usual Dirac algebra; furthermore we work 
with a Minkowskian metric. Note that this choice leads to some differences 
with respect to the expressions obtained with the notations used in [2]. Of 
course the quark propagator is diagonal in the flavour and colour indices and 
the antiperiodic boundary conditions entail: 

S a p{r - (3, x; t', x') = -S a/3 (r, x- r', x') . (A.3) 

Thus the Matsubara frequencies are given by u n = (2n + l)nT. 

In terms of the free fermionic spectral function 

Pf(p) = e(Po)(V + m)5(p 2 - m 2 ) , (A.4) 
the following representation holds for the free thermal propagator: 

S F {tuj n ,p) = - dp - . (A.5) 

J-oo tUJ n - Po 

It may be useful, in some computations, to use the following expression for 
the fermionic propagator in the mixed representation: 

/+oo 
dp e- poT p F (p ,p)(l-h(p )). (A.6) 
-oo 

It is easily verified that Eq. (A.6) leads to the proper spectral representation 
(A.5). In fact: 

S F (iu n ,p) = / dre^ T S F (r,p) = / dre l ^ T / dp e~™ T p F (p){l -n(p )) 

JO JO J-oo 



[ +OO dp p F (p ,p)(l-h(p )) fdre^-^ = - [ + °° dp PF{P0,P) • (A.7) 

J-oo Jo J-oo lU n — Po 



B HTL resummed quark propagator 

In this appendix, following [11], we quote the spectral representation for the 
HTL resummed quark propagator which is suitable for performing loop calcu- 
lations and also helpful in understanding the physical meaning of the different 
terms one gets in the final result. 

In the massless case the quark propagator in HTL approximation can be writ- 
ten as follows [11]: 

*S(u>,p) = *A + (^P) 7 ° ~ 7 '^ + -A-(^P) 7 ° + 7 '^ (B-l) 
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being 



*A ± ( w ,p) 



to: 



-1 

\^-) ln^±2 



(B.2) 



with the quark thermal mass m q (T) = (gT)/VE. Likewise the following repre- 
sentation 



k S(iu n ,p) 



1° — 1 ' P f + °° j P+(w,p) 7 +7-p/ >+o ° p-(uj,p) 



-L 



oo IUJ„ — UJ 



duj- 

oo VjJ„ — UJ 



(B.3) 

is easily proved to hold as well. In (B.3) p± and *A± are related as follows: 



*A±(z,p) = -f du 

J-oo Z — UJ 



00 , P±(uj,p) 1 



p±(uj,p) = -Im*A±(w + ir],p) , (B.4) 

7T 



z spanning the whole complex plane. 

Lumping together the above results one gets the following compact expression 
for the HTL quark spectral function 

Putl(co, P ) = 7 ° ~^' P p + (u,p) + 7 ° + J'V kp) , (B.5) 



with the corresponding spectral representation: 



k S(iu n ,p) = - J 



+ °° , Phtl(^,p) 
du— 

00 IU)„ — UJ 



(B.6) 



for the thermal fermionic Matsubara propagator. This is conveniently used in 
performing loop calculations. 



C Evaluation of G^ s LA (r, p = 0) 



Here we explicitly write the NLA pseudoscalar meson propagator at zero spa- 
tial momentum G^ LA (r,p = 0) split into the different contributions which 
should be numerically evaluated. Recalling Eq. (51) we divide the integra- 
tion over k into three regions corresponding to soft (a), intermediate (b) and 
hard (c) momenta. Starting from the soft momenta sector we have (as usual) 
three contributions of different kind: pole-pole, pole-cut and cut-cut. Their 
expressions are (writing uj± for uj±{k) and f for t/0): 
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2N r r A i 



7T 2 JO 



k dk 



2 (uj 2 -k 



2\ 2 



1 + e~P"+ 
2(3u + f + e -2(3u+(l - f 



2m 2 



+ 



+2 



1 



1 + e-#"+ / V 1 + e+^' 



1 



- k 2 
v 2m 2 



V - A: 2 
, 2m 2 



- ^-) r _|_ e ~P(w+ - w_)(l - f) 



+ 



1 + e - ^- 



2 /co 2 -k 2 V 



+ 



2m 2 



-2(3uj-f + e -2f3uu-(l - f) 



(C.l) 



2N r 



7T 2 JO 



k dk duj 



1 



1 



(,i 2 - k 2 



m„ 



1 + e-P"+ 1 + e _ #" 2 
/5(cj + + cj 2 )t + e -/9(cJ + + w 2 )(l - r) 



+ 



+ 



1 



1 



to 



A: 2 



1 + e - ^- 1 + e~ 
e -/3(w_ + ^ 2 )r + e -/3(w_ + w 2 )(l - f) 



(C.2) 



and 



G (a)( T ) 



7T Z 



Ai /■fc rk 

k 2 dk / / dtU2 — 

J —k J—uji 1 



j-ySwi ]_ _|_ g-/3w 2 



.^(wi, fc)/3+(w 2 , fc) + fc)0-K fc)] 

+ u 2 )f + e -/3(wi + w 2 )(l - t 1 ) 



(C.3) 



The intermediate momentum sector is more involved. We have: 
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2N C /-A2 i2 



7T' JA X 



1 



1 + e _ #"+ 



cj 2 — fc 



•2\ 2 



+ e -2/3w+(l-f) 



+ 



+2 



1 



1 



2m\ 



1 + e-P"+ J V 1 + e+^- 
e -/3(w+ - w_)f + e -/3(w+ - w_)(l - f 

2/3w_f + e -2/3w_(l-f) 
1 1 



+ 



V -fe 2 \ 2 



( uj 2 _-k 2 
{ 2m\ 
1 

1 + e-^- 
cos 4 (a(/c)) + 



2Ar c /" A2 ,2 „ 
+ — - / k 2 dk 

7T Z iAi 

cu 2 - A; 2 



2m 2 



1 + e -/3e fe 1 _|_ e -/3w+ 

/3(e fc + w+)f + e -P(e k + - f) 



+ 



1 



1 



ul-k 2 



2m 2 , 



1 + e-P e * 1 + 
. e -/3(e fc -w_)f + e -/3(e fc -w_)(l-f) | 2 

2N r r A 2 



sin 2 (a(k)) cos 2 (a(k)) 



+ 



7T^ JAi 



fc 2 dA; 



1 + e- fe 



2/3e fc f + e -2/3e fc (l-f) 



sin 4 (a:(A;)) , 
(C.4) 



g Sm^/ a >*/>{tt7 



1 



cui - A; 2 



/3^+ 1 _|_ e -/3w 2 

e -/3(w+ + ^ 2 )r + e -/3(w+ + w 2 )(l - r) 
cj 2 - A; 2 



+ 



1 + e~^- 1 + e~^ 2 



/5-(w 2 ,A;) 



+ u 2 )f + e -/9(w- + w 2 )(l - r) 



+ 



e 

2iV„ Ma 



cos 4 (a(A;)) + 



7T Z 



/■A 2 /-fc X 1 

/ k 2 dk dcu 2 - 3— • • /3+(cu 2 , fc) 



-P(e k + u 2 )f + e -/3(e* + w 2 )(l - r) 



2cos 2 (a(A;))sin 2 (a(A;)) 
(C.5) 



and 
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2N n 



TT* JAi 



A2 rk rk 

1 ,11. I ,1 . / 



k dk du\ 



1 



1 



1 + e-^ 1 1 + e-^ 2 



k)P+(u 2 , k) + p.(u u k)P_(cu 2 , k)} 
-fiiu! + io 2 )f + e -P{ui + ^ 2 )(1 - r) 



cos 4 (a(A;)) , (C.6) 



where a(fc) has been defined in Eq. (52) and e& = \Jk 2 + m^, being moo the 
NLA asymptotic mass. 

Finally, the contribution coming from the hard momenta region reads: 



2iV,. 



7T^ JA 2 



1 



1 + e-P e * 



2/3e fc f + e -2/3e fc (l-f) 



(C.7) 



Clearly one recovers the HTL approximation for G(r,p = 0) by simply ex- 
tending the domain of integration in Eqs. (C.1,C.2,C3) up to infinity and 
neglecting the other terms. 



D Fermion-antifermion loop 



For the benefit of the reader we report some formulas useful for readily per- 
forming the sum over Matsubara frequencies inside a fermionic loop. Defining: 



A(iu; B , E p ) = = £ A s (icv n , E p ) = £ -i- * , (D.l) 

W n + -^p s =±l s =±l ltj V lUJ n - S£j p 

uj n being odd, the following relationships hold [2]: 



^v-a /■ c ■ em ~ SlS2 1 ~H s i E i) ~n{s 2 E 2 ) 

(D.2) 

T ? wAlK ' = W 2 iut - s lEl - s 2 E 2 ' 

(D.3) 

where lo\ is even and n(uj) is a Fermi distribution. We remind the useful 
identity: 

h(-uj) = I - h(uj) . (D.4) 
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